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Abstract 

We analyze low-energy hyperon-nucleon scattering using an effective field 
theory in next-to-leading order. By fitting experimental cross sections for 
laboratory hyperon momenta below 200 MeV/c and using information from 
the hypertriton we determine twelve contact-interaction coefficients. Based 
on these we discuss the low-density expansion of hyperon mass shifts in the 
nuclear medium. 
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I. INTRODUCTION 



The binding energy of hyperons in nuclear matter plays an important role in hypernuclei 
and in the equation of state of neutron stars. The most commonly used approach is that one 
starts from a two-body hyperon-nucleon potential based upon one-boson exchange (OBE) 
and SU(3)-flavor symmetry and then computes the binding of hyperons in matter in the 
Brueckner approximation. It has been shown that uncertainties in the two-body input may 
lead to large differences in the resulting hyperon potentials and in particular in the binding 
of the sigma hyperons |l[]. Therefore it is appropriate to see whether other approaches, 
e.g. a low-density expansion based on vacuum scattering amplitudes can provide additional 
insight. 

Alternatively, Savage and Wise || have analyzed hyperon mass shifts in nuclear matter 
using chiral perturbation theory by expanding in the Fermi momentum. The interactions 
were determined by a chiral Lagrangian with SU(3) x SU(3) symmetry; the hyperon mass 
shifts could be expressed in terms of six coefficients related to the strength of the four-baryon 
contact interactions and the pseudoscalar-meson couplings to baryons and were computed 
at the one-loop level. However, the unknown values of these six parameters prevented them 
from obtaining numerical results for the energy shift at zero momentum. 

More recently, applications of effective field theory (EFT) in nuclear physics have received 
a renewed interest (for a review see So far, most investigations have been devoted 

to the two-nucleon system 0-0], which is characterized by anomalously large scattering 
lengths. The existence of quasi-bound states introduces an additional scale in the problem, 
complicating the power counting required in EFT. A remedy was recently proposed || in the 
form of a selective resummation. The leading-order amplitude of this approach reproduces 
the effective-range expansion with a vanishing range parameter ||. Inclusion of the next- 
to-leading order terms produces a non-zero effective range. 

It is tempting to study whether the approach of Ref. f| can be applied to the hyperon- 
nucleon sector which involves coupled channels, i.e. the scattering length and the effective 
range are represented by matrices. Since our region of application (center-of-mass momenta 
below 100 MeV/c) is below the pion cut, and in view of the problems encountered with 
perturbative pions ||, we use an approach without explicit mesons. Thus the two-baryon 
interaction to leading order involves 6 constant contact terms, since SU(3)xSU(3) can be 
decomposed into the sum of 6 irreducible representations of SU(3). For the s- waves we 
consider, the next-to-leading (NLO) order introduces 6 additional SU(3)-symmetry respect- 
ing coefficients for the terms of order p 2 . The absence of experimental information on HiV 
scattering means that the low-energy hyperon-nucleon data in the present approach depend 
only on 10 (out of these 12) parameters. 

As discussed in the following sections, breaking of SU(3) symmetry by meson masses 
has significant consequences. We model this symmetry breaking by incorporating terms of 
order p 2 coming from one-pion exchange. The two coefficients of these terms we treat as 
free parameters, resulting in a total of 12 parameters. 

In Section II we formulate the coupled-channels formalism, based on effective field theory, 
for low-energy hyperon-nucleon scattering. In Section III we present results of the fit to YN 
cross sections and AN scattering lengths. Some information on the latter is provided by the 



calculation of the hypertriton binding energy [10]. Results for the hyperon mass shifts in 
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low-density nuclear medium and their comparison to results obtained by other methods are 
presented in Section IV. 



II. HYPERON-NUCLEON SCATTERING 

We consider hyperon-nucleon scattering at low energy and use the approach of effective 
field theory. In general the tree-level amplitude in next-to-leading order is written as 

Ao = -C - C 2P \ (I) 

where Co and C2 denote matrices whose elements are the relevant coefficients in the inter- 
action Lagrangian. 

The Kaplan-Savage- Wise (KSW) resummation [KJ gives, in leading order, a scattering 
amplitude A given by 

A- 1 = A 1 -^(fi + i p), (2) 

with M r the reduced baryon mass (for which we take a common value in all channels) and 
/j, the (arbitrary) subtraction point. Using the relation between the if- matrix and the full 
amplitude A, viz. 

pK- 1 = ^-A- 1 + ip (3) 



we obtain for the KSW resummation 

pK-i^-p-^Co + Ctf 1 )- 1 . (4) 
By expanding the term (Co+C^p 2 ) -1 in powers of p 2 we recover the effective-range expansion 

pK' 1 = — + \rp 2 + ■■■, (5) 
where the scattering-length matrix a and effective-range matrix r are given by 

To accomodate possibly large values of the scattering length in the momentum expansion of 
the amplitude one should keep all powers of ap. This leads to the following expression for 
the amplitude in NLO: 



A= -M-„\> i + ip+ M- r C ° 



This is a generalization of the single-channel expression, cf. Eq. (2.18) in Ref. |J, to coupled 
channels. 
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Next we turn to the relationship between elements of the tree-level amplitude matrix 
Co, corresponding to the two-baryon scattering amplitudes of interest. For the spin-1/2 
baryon octet it is convenient to express the 6 real parameters determining the constant 
four-baryon contact interactions through coefficients corresponding to the 6 irreducible rep- 
resentations of the direct product SU(3) xSU(3), which is appropriate for baryon-baryon 
scattering. We write the coefficients of the terms in the Lagrangian for B1B2 ^6364 in the 
form — &s(BiE>2 — > B3B4)// 2 , where / ~ 132 MeV can be identified with the pion decay 
constant and S = or 1 denotes the total spin. Introducing the notation s , Si, s 2 ,t 1 ,t 2 ,t 3 
for the coefficients of the representations [1], [27], [8 S ], [TO], [10], [8 a ], respectively, we list 
the bs values of interest (we note that some of them are given in Ref. [|1J]]) below: 

(i) For the single-channel case S + p — > E + p: 

b (Z + p -> £+p) = 8l , b 1 (Z + p^Z + p)=t 2 . (9) 

(ii) For coupled channels S + n, and Ap: 

b (E+n -> E+n) = \{2s x + 3s 2 ), h{Y?p -> Y?p) = ±-(7 Sl + 3s a ), 
5 10 

b (Ap - Ap) = ^(9 Sl + s 2 ), b (Z+n - £°p) = - s 2 ), 

b (Z + n - Ap) = -^( Sl - s 2 ), C (S°p - Ap) = ^( Sl - s 2 ), 

&i(£+n -> S+n) = |(t a + t 2 + t 3 ), h(T, p -> Y?p) = \{h + At 2 + t 3 ), 
3 o 

1 /2 

6i(Ap - Ap) = -(ti + t 3 ), - = -^-(t! - 2t 2 + 1 3 ), 

2 o 

6 1 (S + n^Ap) = ^=(t 1 -t 3 ), 6i(S°p^Ap) = ^(-ti+t 3 ). (10) 

We mention that coupling between above channels has been neglected in Ref. 0. 

(iii) Similarly there is coupling between channels containing S°n, and An, with coeffi- 
cients equal to those given in Eq. (|lCj), except for an overall sign change for bs{^p — > An) 
and 6s(E°n — > An), compared to bs(T, + n — > Ap) and &g(£°p — > Ap). All coefficients 65 are 
subtraction-point dependent. 

The coefficients of the p 2 term C 2 in Eq. (Q), which we denote by s[, s 2 ,t[,t 2 ,t' 3 (apart 
from the common factor l// 4 to make these coefficients dimensionless) , obey the same 
relationships as in Eq. (|i~0|), since the latter are valid for the full momentum-dependent s- 
wave amplitudes |TT| . Hence they are valid for the momentum expansion terms and also for 



tree-level amplitudes, since the loop expansion corresponds to an expansion in powers of h. 
Let us briefly discuss the breaking of the SU(3)-flavor symmetry (keeping SU(2)-isospin 



symmetry). In a study of its significance it is concluded (TTJ that, apart from consequences of 
the pseudoscalar-meson mass differences and the baryon mass differences the SU(3) symme- 
try is satisfied quite well. The differences in baryon masses we take explicitly into account. 
The symmetry breaking due to different meson masses is modeled by considering one-pion 
exchange (i.e. ignoring the kaon and the eta). To be consistent with the above NLO approach 
without mesons, we expand the exchange amplitude and keep only the leading, order p 2 , 
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term. This basically leads to addition of a symmetry-breaking term —C 2 p 2 / f 4 to tree-level 
amplitude, Eq. ([I]), depending on two parameters, u\ and u 2 . The relevant coefficients in 
the spin-singlet channel are the following: 

C^(E+n -> E+n) = 2u u C' 2 (E + n -> E» = -2>/2ui, 
C^(E + n -> Ap) = 2^/2/3^, C£(£°p - Ap) = A^. (n) 

The spin-triplet coefficients are obtained by multiplying the above values with —1/3. For 
the channels S~p, E°n, and An the only difference is in the sign of the term C' 2 {TPn —+ An) 
compared to C 2 (E°p — > Ap). 



III. APPLICATION TO HYPERON-NUCLEON SCATTERING DATA 

We now turn to the available low-energy data on hyperon-nucleon scattering. We consider 
only scattering with laboratory momenta below 200 MeV/c, i.e. center-of-mass momenta less 
than 100 MeV/c. There is only a small amount of relevant data, corresponding to total cross 
sections for E + p — ► E + p, Ap — > Ap, E~p — ► E~p, E^p — > E°n and E~p — > An scattering 
(see, for example Ref. |[2|). ^ nas been known for a long time [|HJ that these available data 
do not allow a unique effective-range analysis. The most precisely determined quantity is 
the "capture ratio at rest", 

1 3 

r R =^r + -r li (12) 

with 

= gg(£-p -> E°n) 

S as{Z-p^Z n)+as{Z-p^An) , 1 J 

where the cross sections are taken at zero momentum and S denotes the total spin. 

When fitting the calculated cross sections to the data we make two corrections, following 



Ref. [11 1 . First, the S p — » An transitions are kinematically suppressed relative to E p 



E°n, as a consequence of the large momentum release, which at threshold is 290 MeV/c. 



According to |yj this necessitates inclusion of a correction factor r> ~ 0.16 in the calculated 
cross section for — > An. The second correction concerns the non-negligible contribution 
of the 3 Si ^ 3 D! transition in the experimentally measured cross section o"i(E _ j9 — > An). 
Defining the ratio of 3 Si — > 3 D! to 3 Si — > 3 Si cross sections for E~p — > An, a value of 



w 0.57 was estimated in Ref. 11 



We first attempted a fit of 28 data points (the cross sections shown in Fig. 1. and 
the capture ratio tr = 0.468 ± 0.010) with the leading-order calculation involving the 5 
parameters si, S2, h, t 2 , t 3 , varying the value of the subtraction point \i. In this leading- 
order calculation we could not achieve a reasonable fit. 

Then we turned to the next-to- leading-order computation, which added the 5 SU(3)- 
symmetry respecting parameters Sj_, s' 2 , t' 2 , t' 3 , as well as the 2 symmetry breaking ones, 



U\ and u 2 . As further constraint we used the existence of the hypertriton bound state [10 



which is compatible with the OBE potential NSC97f of Ref. ||12|| . The advantage of this 
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p~3|l is that while the scattering data are mainly sensitive to the 3 Si AiV interaction, the 



hypertriton is more sensitive to the 1 So AiV interaction. Thus we imposed as requirement 
on the fit that it had to lead to A-nucleon scattering lengths of a^ N ~ —2.5 fm for the 
singlet, and a^jv ~ —1-75 fm for the triplet state, respectively. 

Several fits of the same quality, x 2 being around 13 (for 16 degrees of freedom) could be 
obtained. If a fit was obtained for a certain value of the subtraction point fi, it is possible 
to vary \i in the range of 0.1 GeV to 0.3 GeV, leading to practically unchanged scattering 
amplitudes and thus the same medium effects. We do not consider members of such a 
family of fits as different. The slight change in x 2 corresponds to variation (in a certain 
range) of the parameters of terms proportional to p 2 in the tree amplitudes, thus changing 
the effective-range values. This is not surprising considering the rather large error bars on 
the measured cross sections. As a consequence, the parameters s[, s 2 , t[, t' 2 , f 3 ,ui, u 2 are not 
very well determined allowing fits with different effective-range values. Our experience shows 
that all these fits lead to quite similar scattering amplitudes and predictions for in-medium 
mass shifts. 

Our strategy in dealing with these fits is to accept only those which do not imply spurious 
(quasi-)bound states ( "hyperdeuterons" ) near threshold. For example, a fit with the lowest 
X 2 led to negative effective range in the spin-singlet Ti + p channel, = —2.32 fm. Together 

with the corresponding scattering length a^+p = 2-64 fm it would imply (14| a bound state 
with binding energy of 3.3 MeV, which is known not to exist. Therefore we rejected this fit. 
With x 2 — 13.5 we obtained a fit with reasonable values of the range parameters and 



AA" scattering lengths close to values of model NSC97f in Ref. |12|]. The corresponding set 
(denoted by A) of parameters is the following: \i = 0.15 GeV, s\ = 1.00, s 2 = 1.78, ti = 
1.21, ta = 3.43, t 3 = 1.87,5; = -0.447, s' 2 = 1.82, = -2.42, t' 2 = 1.05, t' 3 = 3.60, «i = 
—0.260, and u 2 = —0.0998. The resulting cross sections are shown in Fig. 1. The calculated 
value of the capture ratio is r# = 0.4670, the AA^ scattering lengths and range parameters 
are: a^ N = —2.50 fm, rj^v = 1-61 fm, a^ N = —1.78 fm, r^jy = 1-42 fm. For S + p — > H + p 
scattering we obtain then: a^l p = 0.55 fm, r^°) p = 0.36 fm, a^+ p = 0.94 fm, r^l p = 0.35 fm. 
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These values are quite different from those of model NSC97f. 
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Fig. 1: Results of the fit (line) to YN cross sections with parameter set A. 

Our results high-light one of the main problems with the available hyperon-nucleon 
scattering data (apart from the large error bars): the absence of truly low-energy cross 
sections. Fig. 1 shows that the lowest hyperon (laboratory) momentum is larger than 100 
MeV/c, which means that the reciprocal of the scattering length and the range term, ^rp 2 , 
can be of the same order. This leads to results for the scattering lengths and effective ranges 
that are not unique. 

We mention that we could not obtain a satisfactory fit to the data with both A-nucleon 
and E + -proton scattering lengths being close to the values given in Ref. |T2 . 
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IV. HYPERON MASS SHIFTS IN THE NUCLEAR MEDIUM 



Since we have dealt only with the low-momentum expansion of the hyperon-nucleon 
interaction, we can only consider hyperon self-energies in the nuclear medium in a low- 
density expansion [15|]. A simple expression relates the in-medium mass shift (i.e. the 
self-energy at zero momentum) to the hyperon-nucleon scattering lengths H, viz. 



AMy = — 

M r 



i {0) 

l Yn 



3a (1) 
oa Yn 



Pn + 



l {0) 

l Yp 



+ 



3a 



(i)' 

Yp 



Pp 



(14) 



where M r is the reduced mass and p n (p p ) is the neutron (proton) density. 

A better approximation is obtained when one takes into account the Fermi motion of 
the nucleons, which amounts to integrating the momentum-dependent forward-scattering 
amplitude in the nucleon Fermi sea. In this case one has 



AM- 



Y 



2 [ PFn (4°i(0,p;0,p) + 3A? n (0, p;0,p) 



+ 2 [ PFp (A?l(0,p;0,p) + 3A? p (0,p;0,p) 



d 3 p 
(2^)3 

d 3 p 



(27T)* 



(15) 



with (A^) being the singlet (triplet) amplitude, p is the nucleon momentum and that 
of the hyperon is zero. 

Since our fits to cross sections extend only up to momenta of 200 MeV/c, that leads to 
a maximum density of 0.4po for the case of isospin-symmetric matter and 0.2po for neutron 
matter (po = 0.17 fm -3 is the saturation density). Results for hyperon mass shifts for the 
parameter set A are shown in Fig. 2 for isospin-symmetric matter (a), and pure neutron 
matter (b). 
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Fig. 2: Hyperon mass-shifts in isospin-symmetric matter (a) and neutron matter (b), 

calculated with parameter set A. 



The results for the A mass shift in Fig. 2(a) are similar to those obtained in Refs. [JIH-IS 



using the Brueckner-Hartree-Fock (BHF) method with the OBE potential from Ref. |20 



except that there the linear behavior holds up to 0.2po, after which nonlinear behavior, 



S 



leading eventually to saturation, sets in. For the E mass shift the results differ, since in 
Ref. [16] a small mass decrease is found at small nucleon density. We mention that the 



hypertriton calculation [10| shows that the hyperon in is a A with probability greater 
than 95%, thus the behavior of the E in nuclear medium has a very small effect on the 
binding energy. 

Results for neutron matter cannot be compared directly with those of the BHF method 
(and potentials used), since only larger densities were considered in Refs. JIB]-|I%||, as well 
as in Ref. |21|. Low-density results [22] show linear behavior until about 0.2p and striking 
differences for different potentials. Mass shifts based on parameter set A, Fig. 2(b), show 
similar behavior to these BHF results [IE - 15 , 21] , with the difference that the latter predict 
the E° curve below that of E~. We mention that in this case, because of isospin-symmetry 
violation in the neutron matter, the nondiagonal mass shift AMas is nonzero and quite 
large for our parameter values, leading to strong repulsion of the A and E° curves. Without 
this repulsion we would recover the "natural" ordering (in increasing mass shift) of the 
E + ,E°,E~ curves. Also, the BHF calculations seem to prefer a mass decrease for the E- 
hyperons (except for E~ and, depending on the potential, also for E°), compared to our 
repulsion. 

We conclude that due to the scarce scattering data the predicted hyperon mass shifts in 
the nuclear medium are somewhat uncertain, but the uncertainty can be partly eliminated 
using conclusions from the existence of the hypertriton bound state and imposing a con- 
straint of reasonable effective-range values, thus avoiding spurious hyperon-nucleon bound 
states. Computations in the BHF scheme provide more reliable results for larger densities, 
but they suffer from serious ambiguities stemming from the use of different potentials, all of 
which describe the scattering data equally well. 
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